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Anomalously high and sharp peaks in the conductance of intrinsic Josephson junctions in 
Bi2Sr2CaCu208+i (Bi2212) mesas have been universally interpreted as superconducting energy 
gaps, but here we show they are a result of heating. This interpretation follows from a direct 
comparison to the equilibrium gap, A, measured in break junctions on similar Bi2212 crystals. As 
the dissipated power increases with a greater number of junctions in the mesa, the conductance 
peak abruptly sharpens and its voltage decreases to well below 2A. This sharpening, found in our 
experimental data, defies conventional intuition of heating effects on tunneling spectra, but it can 
be understood as an instability into a nonequilibrium two-phase coexistent state. The measured 
peak positions occur accurately within the voltage range that an S-shaped backbending is found in 
the calculated current-voltage curves for spatially uniform self-heating and that S-shape implies the 
potential for the uniform state to be unstable. 

PACS numbers: 74.50.+r, 74.25.Jb, 74.72. Hs 



I. INTRODUCTION 

Intrinsic Josephson junctions (IJJs) in the crystal 
structure of the high-temperature superconductor (HTS) 
E^S^CaC^Os+a (Bi2212) exist along the c-axis be- 
tween each pair of neighboring Cu02 bilayers, with the 
Bi2Sr2 04 layers acting as the insulating tunnel bar- 
rier. Kleiner et ali described the novel features of such 
Josephson junction stacks and continuing excitement 
about research prospects has led to an exhaustive litera- 
ture on the current-voltage characteristics, I(V), of IJJs, 
especially on sculpted mesas, 2 . 3 . 4 . 5 . 6 . 7 . 8 i 9 i 10 . u i 12 i 13 i 14 i 15 i 16 
While it has long been recognized that self-heating 
of such mesas is a potential problem, the observa- 
tion of high, sharp conductance peaks has been in- 
terpreted as evidence for minimal thermal broadening, 
and the peaks were thus ascribed to the quasiparticle 
gap feature expected from a stack of superconductor- 
insulator-superconductor (SIS) tunnel junctions. Impor- 
tant conclusions about the magnitude and temperature- 
dependence of the superconducting gap, A(T), along 
with inferences about an extrinsic pseudogap in Bi2212 
have been drawn based on IJJs exhibiting these sharp 
peaks&^ i 10 : 11 : 12 It has also been argued that the sharp 
peaks are evidence that the c-axis tunneling process is 
strictly coherent (see Yamada et al4). 

Here we demonstrate unambiguously that similar 
sharp conductance peaks are not a measure of the su- 
perconducting gap, but they are rather a consequence 
of strong self-heating. This counterintuitive conclusion 
is based on: independent measurements of the tunneling 
spectra (and A) on the same (or similar) crystals; a sys- 
tematic control of heating via mesa height (number of 
junctions, N, at constant mesa area); and realistic heat- 
ing models. We find that only the shortest mesa data 
look similar to the equilibrium tunneling I(V), albeit 



with a voltage-dependent mesa temperature, T m , with 
T m < T c at the gap voltage. Intermediate size mesas 
reach T c before the gap voltage, due to their greater dis- 
sipation, and the transition to the normal state causes a 
sharp rise in their /(V) without the conventionally-held 
sign of heating, i.e., backbending in I(V). We explain 
the lack of backbending by instability into a nonequilib- 
rium two-phase coexistence for the transition to the nor- 
mal state. A realistic uniform-temperature self-heating 
model shows backbending of the I(V) is not eliminated 
by including the anomalously large quasiparticle scatter- 
ing rate, r(T), that is taken from scanning tunneling 
spectroscopy (STS) data for Bi2212jil Although uniform 
self-heating cannot therefore explain our data, the exper- 
imental peaks occur in the voltage range of backbend- 
ing in the calculated I(V) where instability might be 
expected, e.g., into two-phase coexistence (normal and 
superconducting) with a non-uniform temperature pro- 
file across the mesa area. 

Earlier studies^iLSiSii^ made reasonable conclusions 
about HTS properties based on data only from Bi2212 
mesas. Such studies often reported /(V) similar to 
Fig. 1(c) which displays a superconducting gap-like fea- 
ture. However, a different story emerges when one 
compares these data to mechanical-contact tunneling 
(MCT)^i2^2L22^i24 e .g., the SIS of Fig. 1(d), or 
Heating is a smaller problem for MCT 
and STS as junctions usually exhibit much lower dissipa- 
tion due to: higher specific resistance; smaller areas; and 
improved heat removal as they consist of only one dissi- 
pating junction, rather than a mesa stack. If one com- 
pares published data on Bi2212 mesas to these traditional 
single junction methods, the only mesas that show similar 
behavior, like the spectral dip feature, ,i&I9^o 4 2i 1 22 1 24 1 2£ 1 26 
have greatly minimized heating effects. These latter 
mesas use short stacks (iV~ 10) with either: intercalated 
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FIG. 1: (Color online) Comparison of tunneling data for 3 
mesas and a single MCT SIS junction (d) at 4.2 K, except 
(a) which is at 26 K. Note (a) and (b) use pristine Bi2212 
and are plotted on the left hand axis while (c) and (d) use 
intercalated Bi2212 and are plotted on the right hand axis. 
The dashed lines estimate the normal state resistance, Rn, 
in each case: (a) 12.8 f2; (b) 30 f2; (c) 460 fi; and (c) 15 fcfi, 
The departure from this line at higher voltages, especially 
noticeable in (c), is caused by the energy dependence of the 
normal-state density-of-states and is also seen in the MCT 
SIN data displayed in Fig. 4. Curve (a) is a 300x60 fim 2 THz 
emitting mesa of thickness 1 (im (JV~ 640). Inset: the 100- 
fj,m soft Au wire used to contact the Au film atop the 10x10 
fim 2 mesas in (b), (c) and in Fig. 2. The hook shape is to 
minimize contact force (damage) from the sharpened tip. 

Bi22122^ to reduce the dissipative c-axis quasiparticle 
conductance; ultra-small areasi (<1 fim 2 ) to improve 
heat removal; or 60 ns current pulses^ ' 13 ' 14 ' 15 to mini- 
mize the buildup of mesa temperature, T m . 

Interest in Bi2212 IJJ stacks has been further peaked 
by the recent demonstration of emission of significant 
THz-wavc radiation^ In order to achieve the needed 
THz cavity resonance and high emission power, large 
mesa volumes were used and those are at odds with con- 
ventional trend in Bi2212 mesa research to reduce heat- 
ing, being up to 10 -10 6 times larger in volume. Thus 
the ongoing quest to understand non-equilibrium (NE) 
and self-heating effects in Bi2212 mesas is clearly front- 
and-center. 

A stack of IJJs, in a mesa sculpted on a Bi2212 single 
crystal, offers a unique possibility for spectroscopy and 
exotic device arrays using well-defined, uniform tunnel 
barriers that are stable over the entire temperature range. 
However, the poor thermal conductivity, very large cur- 
rent density and close proximity of multiple neighbor- 
ing junctions in the crystal structure of Bi2212 have led 
to significant concerns about heating. A striking exam- 
ple of self-heating is found in large volume mesas as a 
backbending of the I(V), e.g., Fig. 1(a) and Fig. 1(b), 
occurring at voltages far below the gap voltage, 2A/e. 
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FIG. 2: (Color online) Conductance per junction, Ndl/dV, 
for each intercalated mesa compared to scaled data for an 
intercalated MCT junction (N= 1). Inset: ratio of FWHM 
to the peak voltage, /3, for each, together with those calculated 
for the data of Yurgens et al.— - inverted triangle - and Zhu 
et ali£- filled circles. 



This backbending was show n ' to be self-heating that 
is dominated by the particular temperature dependence 
of the sub-gap quasiparticle c-axis resistivity, p c {T m ), 
rather than A(T m ) as was found previously in low-T c su- 
perconductors and explained by NE quasiparticle injec- 
tion and/or self-heating effects! 31 ' 32 Compared to Fig. l(a 
and b), the voltages of the sharp conductance peaks in 
our intercalated mesas, Fig. 1(c), are much closer to 2A, 
so the effect of self-heating on A will not be negligible. 

A basic question is whether an /(V) curve such as 
Fig. 1(c), which exhibits a gap-like feature, is in fact a 
measurement of A(T). One of the main results of the 
present work is the establishment of straightforward an- 
alytical methods that allow such a question to be an- 
swered. First we establish heating effects arc present 
in Bi2212 mesas by directly comparing their I(V) to 
those of single MCT junctions, using similarly made and 
doped crystals to eliminate uncertainties in sample qual- 
ity or doping level. In Fig. 2, we compare the differen- 
tial conductance per junction, Ndl/dV, for intermediate- 
size mesas (10x10 fim 2 by 12-60 nm height, i.e., N= 
6-30) with the MCT junction of Fig. 1(d), all taken 
on nominally identical intercalated crystals. The differ- 
ences are qualitatively and quantitatively unmistakable. 
Based on these data, we propose a figure-of-merit for 
IJJ stacks that is a quantitative measure of conductance 
peak broadening, j3, given by the ratio of the peak's full 
width at half-maximum (FWHM) to its position, V pea k, 
such that /3=FWHM/Vp ea fc. A summary of /?- values for 
these data and those of others 3 -^ are shown in the inset of 
Fig. 2. It is clear that when (J< 0.15 the high bias spectral 
dip feature (a definitive signature of superconductivity in 
Bi2212) is absent in our data and that of Zhu et alpi We 
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will show below that the observed, sharp gap-like feature 
when P< 0.15 is not a measure of A, but rather marks 
the transition of the mesa into the normal state via two- 
phase coexistence that is due to self-heating. Applying 
this criterion to previously published work, it is evident 
that some mesa studies, which presumed to be measuring 
A(T), were likely observing self-heating effects instead. 

It is counterintuitive that the data for intermediate 
size mesas (N— 12 to 30) become sharper with greater 
heating (i.e., larger N) as this does not reflect the usual 
smearing-out of the junction I(V) due to heating. This 
is especially true here since smearing should be severe 
in view of the anomalously large quasiparticle scattering 
rate, f(T), for Bi2212 near T C XL Thus our conclusion 
that conductance peaks result from heating defies intu- 
ition and its justification is the primary purpose of this 
paper. Our contention is based on a self-heating induced 
transition into the normal state through two-phase co- 
existence. We will also show that the gap-like feature 
shifts to lower voltages with increasing bath temperature, 
and this effect mimics reports in the mesa literature that 
are interpreted as closing of A(T) at T c . However, de- 
spite this feature not being A, its disappearance at T c is 
guaranteed in our self-heating scenario because there can 
then be no transition into the normal state if T is above 
T c . 



II. EXPERIMENT 

Single crystals of Ca-rich Bi2.iSri.4Cai.5Cu20s+a, 
were grown by a floating zone technique. Intercalation of 
HgBi'2 occurred upon heating these crystals in air with 
excess HgBr2 gas at 230 °C for 16 hours and x-ray diffrac- 
tion confirmed the c-axis lattice constant increased from 
1.531 nm to 2.151 nm, as found previously^ The in- 
tercalated crystals exhibited T c ~ 74 K from magneti- 
zation and A~ 24 meV from MCT, indicating they are 
likely overdopedJ^ Intercalation of HgBr2 between the 
BiO layers reduces the specific dissipation at fixed voltage 
by thickening the Bi2Sr2C>4 tunnel barrier to obtain an 
order-of-magnitude decrease in the c-axis conductance. 
We can fine-tune self-heating by of a single control pa- 
rameter, the stack height, since total heating power is 
proportional to N for a constant mesa area. 

Intercalated crystals were cleaved, sputter coated with 
gold and Ar-ion beam etched^ into arrays of 10x10 fim 2 
mesas using photolithography. Our MCT apparatus, de- 
scribed by Ozyuzer et ali24, is also used to contact the 
gold film atop the mesa with a soft, 100 //m- diameter 
gold wire that is bent in a hook-shape to minimize the 
contact force and any potential damage to the mesa (see 
inset of Fig. 1). This wire is sharpened to a diameter of 
5-10 fim at the end touching the mesa and invariably the 
tip contacted a single mesa of the array. The multiple- 
sweep /(V) in Fig. 3 shows the Josephson current and 
the number of quasiparticle branches^ corresponding to 
the number, N, of IJJs in the mesa. 
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FIG. 3: (Color online) Details of curve (c) of Fig. 1 showing all 
quasiparticle branches and the I(V) obtained after subtract- 
ing the initial SIN junction (half-open symbols). Inset: shows 
this subtraction accounts for the slope of the zero-voltage dc 
Josephson branch. 



Single-junction methods were employed for compari- 
son with the IJ J spectra. Both supcrconductor-insulator- 
normal metal junctions (SIN) and SIS break junctions 
were obtained on an intercalated crystal by MCT using 
a much thicker gold wire with a blunt tip. After collect- 
ing SIN data, shown in Fig. 4, a hard contact is used to 
micro-cleave the underlying crystal leaving a chip of in- 
tercalated Bi2212 on the Au tip for subsequent SIS junc- 
tions^ as shown in Fig. 1(d). Since heating is virtually 
eliminated for MCT, these data provide the equilibrium 
properties for interpreting the Bi2212 mesa data. The 
mesa's top-contact gold film forms an SIN junction in sc- 
ries with the IJJ stack: it is seen as a finite resistance 
(slope) for the zero-bias Josephson supercurrent (solid 
symbols of the inset to Fig. 3). This mesa resistance 
decreases with temperature in a manner that is only con- 
sistent with an SIN junction. To correct for this, our 
I(V) from MCT in the SIN configuration is used to sub- 
tract the top-contact SIN voltage from the measured to- 
tal voltage, for every current. The result (open symbols 
in both Fig. 3 and inset) is the I(V) of the IJJ SIS stack 
alone and a numerical derivative generates the dl/dV of 
Fig. 2. The biggest SIN correction yields a 10% lower 
peak voltage, but our qualitative conclusions would be 
unaltered by neglecting this correction. 

The telltale signs of heating in the dl/dV of Fig. 2, 
even in the absence of backbending in I(V), are: (1) 
anomalously high and narrow peaks in the taller mesas; 
(2) significantly reduced peak voltages from the MCT 
value; and (3) weak (iV=6) or absent dip-hump-feature 
(DHF) that is seen above the peak in the MCT data 
(at - 80 and - 120 mV in the MCT data of Fig. 
2). The DHF is a universal feature of all previous 
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FIG. 4: (Color online) The dl/dV from MCT for intercalated 
(solid circles) and pristine (open circles) Bi2212 are very simi- 
lar and for subgap voltages they fit the momentum-averaged^ 
d-wave BCS model^ (solid line). Above A, the energy depen- 
dence of the normal-state density-of-states is clearly seen. 

MCT) 18 ' 19 ' 20 ' 21 ' 22 ' 24 STS> JL26i2&2L2& and angle-resolved 
photoemission2& studies as well as IJJ data^ii^ with re- 
duced heating. 

To validate the generality of our results, we note that 
intercalation of HgBr2 has little effect on the DOS of 
the Cu-0 bilaycrs as revealed by the similarity of the 
SIN MCT data in Fig. 4, with and without intercalation. 
Both curves are fit by a d-wave DOSi^l 



N s (E,k)=Re\ 1 = ] (1) 

\^{E-,rf-A{kf) 

where E is the quasiparticlc energy, A(k)= Aqcos(2<P), 
(j) is polar angle in k-space and r is assumed to arise 
from the quasiparticle scattering rate. The sub-gap con- 
ductance is well modeled by Eq. 1 and the fit, shown 
for the intercalated crystal, gives Aq= 24 meV and T= 
0.55 meV while the oxygen overdoped pristine Bi2212 
crystal fit (not shown) yields Aq= 26 meV and F= 0.6 
meV, and both use the same momentum-averaging pa- 
rameter^ At higher voltages, Eq. 1 is inadequate as it 
misses the actual normal-state conductance and the in- 
dependence of A: the latter produces the DHF— that is 
seen more prominently at negative voltages in Fig. 4. 

The tunneling conductance of an SIS MCT break junc- 
tion in Fig. 2 exhibits both Joscphson (small peak at zero 
bias) and quasiparticle tunneling. Because SIS junctions 
are less sensitive to thermal broadening, their conduc- 
tance peak voltages are an excellent measure of 2 A. Thus 
A= 24 meV from the SIS peaks matches that of the SIN 
junction fit of Fig. 4. The dip feature located at 80 mV 
is more readily visible in the SIS data but can be traced 
consistently^ to the feature at -50 mV in the SIN DOS of 



Fig. 4. For example, using the difference between the dip 
minimum and the peak voltage as a crude estimate, the 
strong coupling boson energy^ gives approximately 25 
meV and 30 meV for the SIN and SIS junctions respec- 
tively, a reasonable consistency. Notably, data well above 
the peaks coincide for all mesas with N— 12-30 in Fig. 2, 
and the lack of DHF implies that all vestiges of super- 
conductivity are gone. Then Ndl/dV at the dip/hump 
for N— 6 is clearly seen to fall below/above that normal 
state value as expected for strong-coupling self-energy ef- 
fects linked to superconductivity, and in agreement with 
the MCT data»^^^^ 

While the conductance peaks in Fig. 2 for mesas with 
iV> 6 are reminiscent of the superconducting coherence 
peak at 2A, that interpretation is incorrect. First, the 
peak is at a voltage that is significantly smaller than 
2A/e, the equilibrium value from MCT. Furthermore, as 
the mesa height, and thus dissipation, increases, the peak 
voltage per junction decreases and that trend is sugges- 
tive of heating. Note that the ratio, j3, of a peak's full- 
width- half-maximum to its voltage (inset of Fig. 2), is a 
small, almost constant value for N— 12 to 30, while it is 
considerably larger for the N= 6 mesa and the MCT data 
that both exhibit a well-defined DHF. The data reveal an 
abrupt change in (3 by nearly a factor of 6 occurring be- 
tween N— 6 and 12. This large change in j3 signals an 
important crossover from near-to-equilibrium supercon- 
ducting properties to severe self-heating. To make this 
result more general, the inset of Fig. 2 includes data from 
Yurgens et al£ and Zhu et alii which are on mesas with 
A^= 10-11. The Zhu data are particularly useful as they 
utilize a different independent control of self-heating via 
the mesa area. Nevertheless, a similar abrupt change of 
(3 versus mesa area is observed in the inset of Fig. 2, and 
it implies that smaller area mesas better utilize the lat- 
eral spread of heat in the underlying crystal (suggested 
in Fenton et alJ^). For j3> 0.3, both the Yurgens and 
Zhu data observe the characteristic DHF but for f3< 0.15 
this feature disappears. 

Note that recent STS data3£ show that the DHF disap- 
pears as T approaches T c , so its absence likely indicates 
that T m > T c after Ndl/dV returns to the normal-state 
value for N— 12-30 in Fig. 2. In addition, the onset of 
normal-state behavior corresponds to the same dissipated 
power, of - 1.05 mW, for N= 12-30. For N= 6, the DHF 
occurs at power levels between ~ 0.5 and 1 mW, while 
for N— 12 a very faint feature is seen at ~ 0.9 mW. 

To summarize this section, we have shown that sharp 
conductance peaks (j3< 0.15) in taller mesas are not 
a measure of A(T) despite their gap-like appearance. 
Rather, their origins are the sharp rise in I(V) that rep- 
resents transitions of the mesa into the normal state, and 
these are found over a range of N. The following section 
presents a detailed analysis of heating as the cause of this 
transition. 
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III. HEATING ANALYSIS 

Ideal BCS superconductors (s-wave) should always ex- 
hibit backbending in their SIS junction I(V) because A 
decreases as the effective T inevitably increases with cur- 
rent (power). However, for weak heating, backbending 
may be absent due to broadening of the I(V) by the 
intrinsic quasiparticlc scattering rate, r, or by a spa- 
tially inhomogeneous A. But if heating is significant for 
2A(T B )/e, where Tb is the bath temperature, back- 
bending is observed.— For Bi2212, several factors could 
reduce the likelihood of observing backbending. The d- 
wave DOS (Eq. 1) does not produce as sharp a jump in 
I(V) when T is zero, but also the large T, that increases 
dramatically as T c is approached, might eliminate back- 
bending even for strong heating and thus reproduce our 
mesa data. We test this latter possibility in Section III. A 
by calculating the /(V) for Bi2212 mesas using experi- 
mental r values in Eq. 1, but we always find backbending 
for N> 12. Note however, that within the voltage range 
for backbending, the S-shaped I(V) allows three differ- 
ent currents, so the situation is potentially unstable to a 
spatial breakup into hot spots (T^ s > T c ) and supercon- 
ducting regions. The resulting two-phase coexistence, de- 
scribed in Section III.B below, can qualitatively explain 
the absence of backbending in our experimental I(V). In 
addition, our sharp conductance peaks, that would signal 
the transition into the normal state through the growth 
of hot spots, occur for all N accurately within the unsta- 
ble, backbending voltage range of our I(V) calculated in 
Section III. A. 



A. Uniform Heating Calculation 

The electrical power dissipated within the tunnel junc- 
tions of a mesa results in increases of the average exci- 
tation energies of the electron and phonon systems. The 
short relaxation times found at temperatures of 50-100 
K allow one to define an effective mesa temperature, T m , 
for these steady-state electron and phonon distribution 
functions. It has been showi*i£ that Newton's law of 
cooling^ 

T m =T B + aP (2) 

represents a good approximation for Bi2212 mesas. Here 
P = IV, V is the voltage across the stack of SIS junc- 
tions, a is the effective thermal resistance and Tb is the 
bath temperature. For spatially uniform heating one 
would always probe a near-to-equilibrium system that 
would be, however, at a current (or power) dependent 
effective temperature given by Eq. 2. Thus our I(V) 
do not represent any constant-temperature equilibrium 
IeqiV, T) found with negligible dissipation (e.g., by MCT 
or STS probes). To emulate our mesa I(V) data, we need 
to self-consistently determine for each I and V an effec- 
tive T m , for which I, V and T m satisfy both Eq. 2 and 
the equilibrium I eq (V,T m ). 



We determine the mesa /( V) for a stack of TV junctions 
by a straightforward solution of these two independent 
relations among the values of 7, V and T m . We generate 
equilibrium I eq (V, T) from standard tunneling theory of 
SIS junctions using the DOS of Eq. 1 and experimental 
values of A{T) and T(T) from STS dataJi These STS 
data arc similar to the MCT example of Fig. 4, but over 
a wide range of T . That studyil used one unintercalated 
Bi2212 crystal (T c ~ 74 K) and found numerous data sets 
that correspond to the observed variation of the local 
properties across the crystal surface. We used three of 
these data sets, labeled according to their low temper- 
ature gaps, A x , A2 and ^3= 22 meV, 23 meV and 25 
mcV, respectively. These gaps and T c are close to those 
reported here, and are thus representative of our interca- 
lated crystals (recall that Fig. 4 shows the similarity of 
pristine and intercalated Bi2212 MCT data). The STS 
data sets, A{T) and r(T), determine the quasiparticle 
DOS from Eq. 1, which is then convoluted with itself to 
generate the equilibrium I eq (V,T), shown in Fig. 5a for 
the STS data set A x , T x . 

To illustrate this method, we consider a fixed current, 
I. Then the allowed mesa voltages are represented by 
two independent, single- valued functions for each fixed 
I: (1) V x {I,T m )= NV eq (I,T m ), where V eq (I,T m ) is the 
inversion of the calculated I eq (V,T m ) and it is gener- 
ally a decreasing function of T m ; and (2) Vi{I,T m )= 
(T m — T B )/aI, which follows from Eq. 2 and is an increas- 
ing function of T m . The intersection of the two functions 
is the self-consistent solution for the mesa voltage, V, and 
temperature, T m , for each value of I. This procedure, re- 
peated for all values of J, leads to the full I(V). Although 
these solutions represent spatially uniform heating, the 
effective T m will be current (or power) dependent. Those 
T m are shown for a few curves in Fig. 5b that illustrates 
our model calculations for N= 12 to 30 using the data 
sets of Pasupathy et al*il (labeled Ti and T 3 , while T 2 
data, not shown, fall between them). 

That all data sets for N> 12 show backbending im- 
plies that strong quasiparticle damping, r(T), near T c 
is insufficient to smear it out. To test the effect of r, 
a fourth curve is shown in Fig. 5b for N— 19 that uses 
r(T)= 2T^,{T) and the degree of backbending is reduced, 
but not eliminated. Although this backbending conflicts 
with our measured I(V), our measured peaks of Fig. 2 
occur at the same voltages as the backbending region 
or the sharp upturn in the calculated I(V), for all N. 
This agreement is demonstrated in Fig. 6 for various TV 
and for data sets A x , A2 and ^3. When the calculated 
curves show backbending, we plot both extrema to give 
a quantitative measure of the region of instability. The 
optimum agreement in Fig. 6 was found by setting a= 
110 K/mW. 

Reasonable constraints on a can be estimated from 
Fig. 2 and Pasupathy et al.— Recall that DHF disap- 
pears as T approaches T c ^ so its absence just above 
the sharp conductance peaks (for N= 12-30) likely indi- 
cates that these mesas are in the normal state at that 
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voltage. These points, for N= 12-30, are all at a power 
of 1.0-1.1 mW. Since the lower bound on the effective 
T m for this normal state, without a DHF, is T c = 74 K, 
the lower bound on a is ~ 70 K/mW since Tg~ 4 K. 
However, the lack of any residual structure in the con- 
ductance data above the sharp peak implies that we may 
better approximate the effective T m for this normal state 
as that for which the gap fluctuations (or pseudogap ef- 
fects) above T c become negligible. From the STS data on 
pristine Bi2211— used in our analysis, this would be ~ 84 
K. But in addition, intercalated Bi2212 presents a much 
higher anisotropy so fluctuations likely persist to higher 
temperatures and a= 110 K/mW is not unrealistic. In 
any case, the accurate iV-dependence of Fig. 6 is a more 
rigorous test of the model than the precise value of a . 

Unfortunately, the backbending observed in the model 
calculations of Fig. 5b is totally absent in the experimen- 
tal curves. Thus the uniform heating model is incom- 
plete, even though it quantitatively captures (Fig. 6) the 
decrease in peak voltage of the mesa data as a function 
of N. This discrepancy can be rectified if one recognizes 
that intrinsic backbending is a potentially unstable situ- 
ation that can lead to two-phase coexistence. 



B. Two-Phase Coexistence Across the Mesa Area 

In order to understand the very narrow peaks with 
their absence of backbending, we postulate that the tran- 
sition to the normal state can be via two-phase coex- 
istence across the mesa area. One part of the mesa 
area could exhibit a smaller-than-average current density 
and power dissipation to remain superconducting (T<T C ) 
while another portion, a hot spot, supports the normal 
state (Th s > T c ) with a compensating larger-than-average 
current (power) density. The instability into an inhomo- 
geneous state follows from the intrinsic S-shaped I(V), 
calculated above for uniform heating, that exhibits three 
different currents for any constant V in the backbending 
region. The Au film atop the mesa and the high trans- 
port anisotropy of Bi2212 assures a uniform V across 
the entire mesa area, even with two-phase coexistence. 
The highest current solution represents the normal state, 
the lowest current solution is the superconducting state 
and the middle one is likely unstable. Once in this non- 
uniform state, further increases in current require a larger 
hot spot, suggesting that it would then propagate across 
the entire mesa area as the current increases. The power 
density needed to maintain the hot spot at T c should be 
weakly dependent on its size, so the increase in current 
would occur at roughly constant voltage, in agreement 
with our sharp conductance peaks (Fig. 2). 

Such a NE driven state of inhomogeneous current with 
a sharp conductance peak has been previously docu- 
mented in k>w-T c superconducting junctions,— and the 
occurrence of hot spots (T/j S > T c ) with a size propor- 
tional to the current has been reported recently in very 
large Bi2212 mesas4£ However, unlike our data, these 
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FIG. 5: (Color online) Numerical simulations of I(V) using 
A(T) and r(T) from reported fits to the experimental STS 
data.— (a) equilibrium I(V) for various fixed temperatures 
T, using Ai, F\ from Pasupathy et al.— (b) non-equilibrium 
/( V) including self-heating for various N and two of the data 
sets Ai, A and A3, 7~3 taken from Pasupathy et al»i£ (calcu- 
lations for A2, 7~2 lie between these). Dashed curve for 7V= 19 
shows effect of doubling the quasiparticle lifetime smearing, 
i.e., r= 27V For two curves the local uniform T m values are 
shown and they correspond to points on the equilibrium I(V). 
This correspondence was used in Kurter et al.— and Ozyuzer 
et al.,— but in those cases only the linear, low-voltage portion 
of the equilibrium I(V) were needed, e.g., Fig. 1(a). See text 
for further details. 



cases exhibit a jump in the I(V) to lower voltage and 
higher current as an initial hot spot nucleates. Rational- 
izing this difference requires a digression into the ther- 
mal boundary conditions for various mesa configurations. 
That an increase in current leads to a larger hot spot is an 
integral part of any two-phase model and is justified by 
the data of Gray et al. and Wang et aL 31 i 40 The relevant 
point is whether the mesa voltage increases or decreases 
as the hot spot expands. 

The data of Zhu et al. r- who varied only the mesa area, 
arc insightful (see Fig. 2 inset). Decreasing the area led 
to less heating, and eventually, for mesa areas < lpm 2 , 
they found close-to-equilibrium properties. This is most 
readily explained by a more effective lateral transport of 
heat in the underlying crystal as the mesa area decreases, 
and this idea was expressed by Fenton et al^ Translat- 
ing this idea into a two-phase model, larger hot spots, be- 
cause of their poorer lateral heat transport, require less 
power density to maintain (T^ s > T c ). The power density 
is (V/t) 2 1 ' PcNc, where t and p C Nc are the mesa height and 
normal-state c-axis resistivity at T c , respectively. Thus 
as the current and hot spot size increase, the required 
voltage decreases since t and p C Nc are constant. This ad- 
mits the possibility of jumps to a larger hot spot, with its 
lower voltage and higher current, since then such a jump 
can be consistent with the electrical bias load line. This 
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scenario is presumably the explanation of the jumps seen 
in the large mesa3 and our Fig. 1(b), while the boundary 
condition for the jumps in thin films is likely different. 31 

For the intercalated mesas reported here, heat re- 
moval through the underlying crystal is more problem- 
atic. The large HgBr2 intercalant molecules should be- 
have as phonon rattlers^ to reduce the thermal conduc- 
tivity (ratios of >30 were found in Nolas et al4i) com- 
pared to pristine mesas. As such, the relatively low ther- 
mal resistance of our 100 /im-diamcter, ~ 0.5 cm long 
Au-wire (~ 2 K/mW) contacting the 40-nm Au film atop 
our mesas (Fig. 1 inset) is a more significant heat sink. If 
dominant, this heat sink would critically alter the heat- 
flow path from a model that considers the underlying 
crystal as the only heat sink. Then the hot spot should 
nucleate farthest from the (random) position of the Au- 
wire contact. As the hot spot grows toward the Au wire, 
the added regions in closer proximity to the Au-wire heat 
sink need a larger power density to maintain T^ s > T c . 
Thus as the hot spot increases in size, both the total cur- 
rent and voltage will increase. Then the differential resis- 
tance is always positive throughout the two-phase region, 
and that is in accord with our mesa data in Fig. 1(c) for 
intercalated Bi2212. This behavior occurs in spite of the 
S-shaped I(V) for that mesa if the heating was homo- 
geneous. Another way to visualize this transition to the 
normal state is to recognize that the Au-wire heat sink 
produces a thermal gradient across the mesa whenever 
power is dissipated. The most remote points will reach 
T c first and that phase boundary will move continuously 
toward the heat sink as the current (power) is increased. 

The previous discussion explains how two-phase coex- 
istence can eliminate backbending in intercalated mesas, 
but can we understand why Fig. 1 (b) with the same mesa 
area and N— 28 behaves so differently? In curve (b) there 
is a jump seen for IRn/N~ 3.9 mV announcing the for- 
mation of a hot spot^S while the backbending part of 
the curve may be effectively explained by the temper- 
ature dependence of P c n{T) in the regions outside the 
hot spot^ The only significant difference is that mesa 
(b) is made on a pristine, unintercalated Bi2212 crystal. 
This yields an order-of-magnitude larger power density 
for a given voltage per junction but more importantly, as 
mentioned above, a considerably higher thermal conduc- 
tivity for the mesa and its underlying crystal than for the 
intercalated mesa (c). Therefore for the pristine Bi2212 
mesa, curve (b), the thermal agenda may be set by the 
underlying crystal rather than the gold-point top contact 
(Fig. 1 inset). 

For the N= 28 mesa of curve (b), the heat transfer 
coefficient is found to be a= 38 K/mW, by using the 
methods of Kurter et alJ£ and the full set of I(V,T). 
The above estimate of the conduction through the thin- 
Au film contact atop this mesa is much smaller (a~ 110 
K/mW). Thus the underlying crystal is the predominant 
source of cooling for this unintercalated mesa. 

The fraction size of mesa comprising the hot spot, /, 
can be reasonably estimated. First, we find the con- 
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FIG. 6: (Color online) Excellent agreement is demonstrated 
between the experimental conductance peak voltage, V pea k, 
(large solid circles) and the gap-like features of the calculated 
I(V) for each of the relevant STS data sets, (Ai, A2 and 
A3) and a thermal resistance, a= 110 K/mW. Examples of 
these calculations are shown in Fig. 5b. Inset: ratio, j3, of 
FWHM to Vpeak for the calculated curves that do not show 
backbending. The dramatic decrease between N= 4 and N= 
6 mimics the experimental one seen in Fig. 2 (inset). 



tributions to the total current, If, from the hot spot's 
normal-state current, fVf / Rn, and the quasiparticle cur- 
rent in the remaining part of the mesa, (1 — f)I qp , where 
I qp (Vf) is the lowest branch of curve (b) of Fig. 1. Thus 
/= (I f - I qp )/(Vf/R N - I„)~ 0.26. Next, we esti- 
mate the power density needed to achieve the normal 
state, Tjv, using the hot spot voltage, Vf, and the heat 
transfer coefficient, cty/J. Here y/J approximates the ra- 
dius dependence of heat transfer through the underlying 
crystal^ since the a determined above was for uniform 
heating of the entire mesa area. This analysis gives /~ 
((T w - T B )R N / aVf) 2 and its evaluation for /= 0.26 im- 
plies a reasonable value of 84 K for the only adjustable 
parameter, Tjy. This is consistent with superconducting 
fluctuations (or a pseudogap) found out to ~ 85 K for 
unintercalated Bi2212JI 

Although providing a justification for the observed 
jump to be a hot spot, this calculation does not rule out 
its occurrence at a lower current in the backbending re- 
gion of Fig. 1(b). From the above analysis a lower current 
would correspond to a smaller hot spot, i.e., a smaller /. 
The heating analysis by Fenton et al^ implies there is 
a minimum size to achieve T> Tjv in the hot spot. To 
find it, we integrate their formula^ for the steady-state 
temperature rise using the parameters of our hot spot 
and get a h-pm. radius for ST = 80 K, whereas the above 
analysis leading to /~ 0.26 implies a radius of 2.9 \xm. 
To expect much better agreement may be unreasonable 
given our imprecise knowledge of the proper values to use 
for the mean thermal conductivity, extrapolation for Rn, 
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FIG. 7: (Color online) The bath temperature, Tb, dependence 
of the I(V) for the N=19 intercalated mesa. The gap-like 
feature at ~ 33 mV for Tb = 4.8 K is seen to decrease with 
increasing Tb and completely disappear for Tb= 80 K. 



Iqp, GtC. 

To summarize, the variety of behavior seen in Bi2212 
mesas (Fig. 1) can be understood by a coherent heat- 
ing model that includes two-phase coexistence for highly 
driven mesas and uniform heating for moderately driven 
mesas. 



IV. ADDITIONAL FEATURES OF THE 
EXPERIMENTAL DATA AND MODELING 

Although incomplete, the uniform heating model ac- 
curately predicts (see Fig. 6) the measured conductance 
peak voltages of Fig. 2 as the voltage of the backbend- 
ing instability for TV> 12. This optimum agreement uses 
a reasonable value of the one free parameter, 110 
K/mW, and thus provides convincing, supporting evi- 
dence that the sharp peaks ((3< 0.15) are a consequence 
of strong self-heating. Since the agreement in Fig. 6 cov- 
ers the range from weak (T m <g; T c ) to strong self-heating, 
it is of interest to see if the calculation shows the abrupt 
change in (3 with TV that is found experimentally. The 
inset of Fig. 6 shows the (3 values derived from calcula- 
tions for small TV values for which backbending has not 
yet developed. Note that the abrupt decrease in the cal- 
culated (3 between TV= 4 and N— 8 is reasonably close 
to the experimental one (Fig. 2 inset). 

This abrupt decrease in (3 signals a crossover in how the 
tunneling conductance achieves its normal-state value: 
for large TV it is due to strong self- heating ((T m > T c )) 
while for small TV, the mesa remains superconducting 
and dl/dV merges with the normal-state conductance, 
ijv, as it does in thermal equilibrium when V equals the 
value of 2A(T m )/e, where T B < T m < T c . This perspec- 
tive comes from comparing the calculated I(V) for N~ 



4 and TV= 8. For TV= 4, T m is only 38 K at the gap volt- 
age, V= 2A(T B )/e= 48 mV, while for TV= 8, T m reaches 
T c = 74 K and the conductance reaches Yjv at a lower 
V= 44 mV. This analysis affirms our previous assertion. 
Thus the model confirms what the data already told us 
and it distinguishes (3 as an appropriate figure-of-merit 
for strong self-heating effects that mimic the supercon- 
ducting gap feature in I(V). 

Our data, modeling and discussion, so far, have con- 
centrated on Tb~ 4.2 K^C T c . Figure 7 shows the de- 
pendence of I(V) on Tb for the mesa with TV= 19, for 
which we have already established that its gap-like fea- 
ture, shown for 4.2 K in Fig. 1(c), is due to two-phase 
coexistence caused by self- heating. Upon increasing Tb, 
this feature shifts to lower voltages and disappears com- 
pletely at T=80 K where a nearly linear I(V) is found 
out to y~ 40 mV. Thus the gap-like feature appears to 
close at T c in the manner of the superconducting gap 
in conventional superconductors in equilibrium. How- 
ever, such behavior is also expected from self-heating. 
As Tb increases, less voltage (heating power) is required 
to reach T c and for Tb> T c the mesa is already in the 
normal state so no transition is possible. This discussion 
points out how misleading the heating effect can be. Not 
only does the sharp upturn in I(V) mimic the expected 
behavior of a tunnel junction at the gap voltage, but its 
disappearance above T c seems to (incorrectly) confirm its 
assignment as a superconducting energy gap. 



V. SUMMARY AND CONCLUSION 

In summary, we find that the occurrence of anoma- 
lously sharp conductance peaks in our data on Bi2212 
mesas is an unexpected, counterintuitive consequence of 
heating. Guided by the equilibrium spectra of single 
MCT junctions, these peaks are explained by a tran- 
sition to the normal state through two-phase (normal 
and superconducting) coexistence, and it occurs at V< 
A(Ts)/e. The data also provide a measure of the sever- 
ity of mesa heating by the figurc-of- merit, (3. Our uni- 
form heating model, using values for A(T) and -T(T) 
from STSjiZ shows quantitative consistency with the TV- 
dependence of the conductance-peak voltages and the 
abrupt change in (3 between TV= 4 and TV= 8. This, to- 
gether with further analysis, gives strong evidence for us 
to conclude, unambiguously, that the sharp conductance 
peaks {(3< 0.15) are not a measure of the superconduct- 
ing gap, but they are rather a consequence of two-phase 
coexistence due to strong self-heating. 

Sharp conductance peaks are often seen in intermedi- 
ate size mesas and they can be easily, but according to 
our analysis, incorrectly, assigned to A(T), even though 
such two-phase coexistence will mimic the closing of the 
gap-like feature at T c . Thus the sharp peaks should not 
be used to make inferences about A(T) near T c . As 
a heating phenomenon, the sharp peak must disappear 
above T c : this behavior is easily misinterpreted as a clos- 
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ing of A(T) and we reiterate the important consequences 
of such a false interpretation. If a measurement seems to 
show that the superconducting gap closes at T c then it 
must necessarily lead to the conclusion that any pseudo- 
gap observed above T c is extrinsic to superconductivity. 
Such a conclusion would be at odds with other spectro- 
scopic measurements . 17 ! 36 ' 38 

The extreme difficulty of eliminating heating may 
imply the need to reinterpret some recent IJJ stud- 
— — — — — that generally exhibit small or non- 
existent DHF in dl/dV and (3 values of 0.03- 
0.15i 8 ' 9 ' 10 i n ' 12 In view of the strong heating effects, one 
may also want to revisit the interpretation of the broad 
peaks seen in mesas near and above T c . A comparison 
of the present work with other IJJ studies indicates that 
when (3> 0.32 heating may not exceed T c in the mesa 
at V= 2A/e. Then the mesas are closer to equilibrium 
and the upturn in I(V) is a measure of Z\, albeit at an 



unknown T m somewhat higher than Tg. In this case the 
voltage dependent T m ensures that /(V) will never per- 
fectly replicate any equilibrium I(V) found by MCT or 
STS. Another indicator, beyond /?, of strong self-heating 
is the absence of a dip/hump feature, a characteristic 
signature of superconductivity in Bi2212. 

While our conclusions reduce the usefulness of Bi2212 
mesas for fundamental studies, it is important for the 
scientific community to recognize this limitation. 
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